“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1993-03 


Two-dimensional boundary surfaces for 
axi-symmetric external transonic flows 


Al-hashel, Waleed Isa 


Monterey, California: Naval Postgraduate School 
http://ndl.handle.net/10945/24243 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sha Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 










































































































































































































































































































































































































































































































































































































. -s6eoe8 Onto Gal cae re =" * « ae” ob ig 4, aaatan * fing? Op 
é 4 Ss r 4 me x ep 8 , t ecelr ies i Pal 
* aan * is e any Geese aa | es seeeep Ca ee Se tes appa se <0 8 lobbies 
a 5 ‘ " jisle-ey qe wal eter Soe Tapia eta ph ern Ai eee me » 918 9 a3 i gesierace 
i =e ote we . Pe] 4 ‘e 4ete 
« : Fs Ree ad ena ie ada Oris is : 1s oto Hee sueeets nyt in é a ratiad 
‘ 
: cron ‘ ' 6 4 
. Fi ee a . ¢ 
oot + a i po 4 ‘ ‘ 
3 aoe ' Tae iy A Je 
Ft . eo ree ‘5 ’ t 
a ' . + ' ’ Tit fete? i bias 
8 : yo : 4 ‘ , oe 7 De ears: “¢f oe ake cH 4 Ve edad eancees a 
2 A ‘ a . ’ eo. 4 S98 46 ahead ee he et eee 
‘ ’ ® hia. i OF 
ne ad ase} “58 nt ice 
- rsh 
* é » . . 
a fs ’ ‘ . 1 oe orig " 7 
ane an) ® é RC fut 
‘ a A 
, ee Mpcaedy 
. rt tas 4 roe aie “gi Hy gestae re 
‘ : ; sat de aA iegt pee alacant aia 
ian oe ses ta 20.8, here doh ¢ age 9. aety 8) RNY on iD 
th Ne ae 2gtas¥ gad Car regres Feta red ees a Was ae rete Pye 
. eee s . ¥ waar ir ay ac ie ¥ hae oy tb Ae ee! ih Vu ed Saba 
ee ‘ ra Oe Ad ty ats pute te rich asee » at air £44 wend a eta ® 
f. 2 sre . “ve aah aaa “oe be a) tn 
‘ s Part » ’ . ic hos £5 ee Fi * as 32> 3 ras oA 
‘ ’ 4 . Coe sh ty A wi. “ine Be eeta te or P iyi 
ak "h8 
‘ . ret ‘ are 3 : hg td ast. 
; ° ’ * ‘¢ *O yd ob 5 450 pte ot veg 
‘ \' . ‘<n oF : ‘ een Bane 9 & 0 Sgd ne kie b 
. * >k : : ay Ue 4 
' . ‘ ’ ate 
: eat te A ¥ spetg ky 
+ tee a ‘ > are D eiveae iad 3 9 ana eta 
ve 4 oe i atghat, a tlarkye. Be Ngreh ae aba vay) re 
i ake Ly 
, : : ny ey: jae ety as Speake tsetse 
: ; : ¥ Se ali fe age sh Aah a5 
wows ‘ . a aa ate et he a me a oa ha heer 
ee Pebts ” Ps ue | Cr ay Wes wae FLA Meat tem It PP Et J Te bs ie beth oe yD Ate trade ay sh Ps! oF) “ay 
. . . ea > Ree, Beste on ok J ee ’ ref ot a gett baeehy Mg 33 + are | rs Wd Re Saad she bax coaneanl 
' soe ¢ Aaah ete Leet F828) c ace Sar? fre" a va we vents =f, oe re oe epee 
en rerio : ae > ae ofa tps Se tas we $6 4.Raa righ, Lae! fag raat: vist sg tbe 
Cetra et Ta “ ' ae P , i F . Busoni : ; SeOTiee Nalitse ackuneceee Fie) eer “¢ Fig LA Ree shies as sat 
. ‘ ’ 
ies a Suess ‘ . 8 yaa ; 1 sod eo ome . 7P oe 4 , re ire eo ved, “Worse 8 hgh, AP ire radia a 
. = . ve . > . <= one tens “aoe Chet. 5 » 1 i. H i‘ . ‘Pad Seb ‘ Cay ry shee sR 
yee ' Ae . Vi ’ ‘ ¢ Pe | ‘ . 4 PU ibeeh is Pie) * we ere eb pity : ae Nas ry 
Se a o> oe » eu D ee ” a? Sooty ate beat ew ore ae 1 “2 & 
1 , , ' ’ ‘ ’ ‘ » CF Se aa Ag ‘ ‘ » * ‘ . "> . 298d nk eee) a* Re 
’ ‘ bo owe) he ON ‘ , te oo . ats ’ F ou a 1. vy : ity 
a ‘ rt ' . ‘ . 8 . > ie > ar ey a sMeo 7 ue 
Mi : > Bae “ os tr , , a wi eage . . +e on id f,s 6 “ ' t ¥ gfe . Ce Pd no uae ® Oya eer. ld wf Vg Ae tab eighty 4 
‘ ue . eo : - we wt . sy 1s at oe , . ,*  ¢ . PS he Te De aT ‘ ata Cc gah heh Bet a ails ‘tA Baw, a os Bi Praseete % fod dal ed, al 
ees uM a at “ . aes i + Fe Pay) ‘ ewan) Pome} ' eM Peli mer recs ore ir) a ® r yh Eola eae? ‘ i Nrb 1, ae 5h. 4 oteuno, Spreng sh tee hiss Shy, + ALN tates 
Sera aR ! P ‘i taeewee wdte H 100 ‘ “ ri Pnicy wees She: | lerereby rk frm) (ei beparg ee eo 4 u 4 waa! > Lan wep eeeh? D ih De bgad er an SM LE WOE sea ey Be ited a6 
i : Pas von ’ ‘ i Seen Tced fo i ' oe se oUt ree wy ORE = ye a) Pere ele beear, sad iss emang ope elgg a be TaAY) As teh 209 bat oy oye 4 
’ as ‘ % a BE eee erat i e mies ’ . «fA bor de ano of Fe oy Tada 7 “ Vee 0 Grd tomas 4 nae LP) aM, aby fig We tk Pa vate WS gb ied ye J eB ou 
. : a i te Came: >wt - . , oo, Wiese a testa ate a” r 
a te + <a e Saws akaioe F aes ht Ms . P Oe ee eer moa Perro ke ce ‘ j Eee a. 
* ae , ve : ’ . , Z vue oo ety . vf E ‘ ‘ ety onde Lays * le 
. * . ae te . . t¢ . or . ' sie = we set et "8 ‘ ‘ oh 
es ass Aas 2 a ee pao eer eat tee PY oarak (nnn Sear Ya Mere ee Re taere ogehe te gion he 
’ ‘ i ' ta ‘ t ie, uty PLE Sy ' ' - , s* " 6 erode ot ef . 
« + @ Mer eee | ’ ‘ an foo . Pe Ye. er | ‘ . ' 2 . ‘ > ae 
; i a 4 ay . * iy oe . ‘og ‘ ete HI Ew nee i cp. ‘ 1 tw - sees 
ay t Ja sQ eageieye 
: a 
ne H idrerdpe ny, A shah 
+ 5 ‘ ar ace a. , as " a % 2 
Ou Thea here UAE daed tier aay a a whee A F att ees 
. 4 ' 3 ee a oan Nias sant med PL Veter 4 
fore ohh wth 1d, etiea Fee 4 ert oe $8 wale teas eri i ae (aia ach het 2S aa Ar 
. bree ah bow ob , ayes aye ugh a ere eee ee ON te ah, Saas shan VM WS Bae be ba bom 
Re aURLaeelaaeaw en iin, Pie Se 348; ag wildersisbas ue sai ae ane bitakob ha here anal 
aud tay a ia Bett FG og bp: aay ase rie Tse * rans hi 
od vtebebeb GS Rye te gynst PUP RNS Tir 425 egies Aa AS = bent aw 
1 ‘3 fertse seed. wih % us + cee aeeeente vr) es, oe ‘~ 
’ ’ ate: Arild eas Peporer 
w 
_ é| 7 fauha ; 
. . Jitte pba aA > tet ‘ey wi dyrs 
‘ phe a5 Hast weedy S4gh 8 jsla* 55 
F V4 calle dy ‘g2d08 aavarer ere aaa ft 
o* z Warrant 9, \ “a Vaya! d ates 
ef vay ag he : tie k ’ i ang het 
. Posty ere tesery Genus seh hae aoe <s ated béoad 2 Vs Wass 
"talat't * @ gt & oir yon eé wa ar a tye mein tay 
¥ ws e9a gto 4 ohod syere PL ae ery z, B Phy teen oe wi Newa ayes eae 
4 Bs 6 eu oe ee ee ae Boiwiree oy iti % tb, Nines 
_ ae | gare raed Lei er 14g Tes Feil “a yy 
" , . . to tes ey 
42 8 , aoe rout, oy ‘er ys 
rs ’ , <a ks ’ . t . bi ‘ + a yfere Boeke > * eas 
‘ : ay t sii ane ' Re ae the eartte "i arene ats £ 1 . Peak ahr rt yee ty sigh” abs ayed ot : 
A g ber - An wes +4 ', me Wg 4s rs » ’ 3 
« en : ¢ Uns nen Z H Ua gas thas " i ny yr oeras an 24 ay Le 1 if ray a rf oy tans » . oy te tee ‘ i , my 
’ ’ : . . te as 1 ‘ ae eleo ‘ . sv thee = tha oy ? o> 6 by 8 aay 
. ‘ o 4 ts ; a PG . ae eae ae rae Peer Sel or ems aa A Ls ? . P Pahl eaacelat ee GRRE RS 
. . sa ® 1 5 aaa Nd “e 84 eee) he art eg a a er Le ae he ¢ ‘ . ated 
' ‘1 wa eon ' r Car » ree 5 mee ‘ Tee ‘ 4G 8 PP set 4 vf tea pe 
. e ogt . . 2% oe ieee i bie 
ry ace A ‘ ‘ . . ¢ Py eh} Fi wd re, Se ‘é ry Peet Ty) amet 7 ' ‘ ’ hey ae TOE vs mee . 
’ ‘ . 7 * a8 bata gs ". eo + ‘ t 2 *%abeet . i * eye 1 a *s » ’ 
. : ie gh ; ana a . Fa Fa , os + An oe uae Wie oe ’ eae Gia nutes *s «  &e'* . ‘ A Ae ne 
a? » ' . = . . m4 - s ‘ ademas “ . oi ip se ntay Gah! 
. 8 » ‘ a a ’ . 
: Fae t : Moers er 5 aaa re ao ‘ »* 1 "oN beagh te 
ong 3 . oe . ‘ "es + ‘ Fy feb 2 Prateek 22 2 | ‘se 
a ? o 
3 os H F . | ; | s*ee 
F . ‘ ’ ah oye 
A vs a * * F) oo. é te * ue’ 
LS « aie ay ci Merhe at ; Stoke thet en ae 
sd ' oe o* 4 Hee i a] . < ‘ fe et 
: ~ i : oat svat + ti ' : t ae ae le ’ > #0 Ot a Oe ‘ 
: Ea " 1 oo ay ote leh 2 8 8 CU e ptt e be } 
. re Bate att “se bs + *eeten I 
F na ren aaa Marre aS rise ae LT ee aT] i pdeyenghel wees 
‘ . * . - 7. % ‘ (ae) 
Cen By E ; ae ; : r riieneries ats 
4 sf tes “ ‘ ‘ i vabas 
re ‘ ’ se * bee - . oe ee \. 
' « r+ hae oe so 4 +e an 
‘ * ft ete * ‘ . ’ ten 
p a - 4 ‘ 2 3 PORCH Ul 4 Dat] ; 
. a sf a! ae td , Pic’ ‘ 
7 one 4 + ae + 
ae in Wyonae Sih aens a} Shee # 
‘ . Sieut . ¢ a ac me: . ri 
oye one " a) is 5 . 
" * 4 6 1 rh. te + 
Celt (Peete peed) has i fat ef . ° « 
Pr + ‘ ‘ ate ‘vt & > rs ee 
. Ler ie ES . Balt rca aie | we Ae a Ys . . as * 2 oe dntys org 
> rae ry « » ’ 4h etre , eh tte ait 1 # whee s Lae Sd yes ate rn the 
. 4 » Py er “ae -.# 1 og » 











¥ 
ee) Le a 





Feet. or 


“ah at oh 
- 
f 


Pae%a2 






« 


i 


Hhotig en” 


art 





















yas al sogtet « 


Ay tes. Ab ws 


ae AS 

ei o fet tee 
enh see 

ake 


hs Oh hat il ” ai 


0% pins 
Aves ay ost. &3 
ie we SEN as 

935 





ae ete 







* 
ay 


a et De 
ary g?h 























































: f 
pa! 
: : 
Al piety abe yg o 
7) + Phebe mee ah 
e 
rang ay 
; , lH TA 
arsed gay, be PAT iy 
’ y oe thd gh Ete Fg avage See a 
. . ta ale ae ‘ 4 
£ebesry. Lee pwbay lout ae i 
ve ry a \s 
. Pesan ne JM s 
te CY Tate WRthe sets s 2 
Ses ageanerty ri 1 Seett me ot Fy ‘ 
3 ae is be 


Pat heh ae 
A 
















: 
Smee ae p 

Lt agaee? 

= Me peethey aie ? on 
"~ 0, AERA 
? 





“Hg ed ohete tee (ae 
eguae gee at 






















EG 




















































































































































































































+t ae sit Pa we) 
et 1 arye SpreMste gslel ie uaaet re hens ape REE 
‘ . OP) Oe OP aut Hh Org gt ’ Sar ance whryrse a" 
vo EAP LONE TE TEs ge of Aen Aref eae ’ at twee 
2a ve fae rf liye. ue A fot, BRM yt ee FRO ? 2 Arh abet ig 
a Sy “08 “ae fT ty te tac s sent i tesa y we 
Pea a bys upeD ety pt aH." 4's en ites « Fepehe # oon 8 wan en ; ple 
’ . ° ed ’ . ey ’ f ‘ Bete Fe 2 8 oat ayeraa ety © at gg trary er aye were 7 tNs t 
¢ > « ? . ‘ fe ses ees pee ep eee ey Prk uh Soe Tsy pa tae ae yee 
1s ae at ' \ , . cl Ld + s+ H $ teeters : * Ae 4 hit Some bk were} f 
. fy t Po en Siew ong rs oe a ee . he > 4 24 aly ener y | : i Ps GCF gel Se emo, 
es , tee ’ @» . eee re ’ CRY te 3 q s : Padre es 
. rn ‘ rs Sed we . 4 ote shy tee game Hye a ale Pee ad Miietsly f Sas eh heres ¢ aaa 
Pv) ’ Ch ry + aeee ao ao MT esac) I al io flee TaNe Seba abet eee 7” whens? Smeran 
"> ‘ 2 ] oy ar ee: ‘ * , Lares FLA haa a ; maT pl, ‘ed 
4 tere UR ' arya beg iets 40 Ae he ay 
ity Mohr rama fy acct eae a, 2 fobs 
’ > g sh weg . tle * bard * 
4 PG. 44¥ Ml Sie, ey ee) eee ‘s Carte 
. 2.8 = gk pita er heat ry pre i "i , meee ye 
5 (meer Pek ae Pe erly or oer ted our e cm Le et Say a é at ,¢°4€% np Caress “88 
coger PR oe Dees Ra) OC) mo aot bia deny Viete vis sae “es? 
eo Pe es tend rer en) ee a ee ry Th Taigte “s i + aniyaly yar rer 
“eo #8 . eo era Ne 4 r 4 gare aes sheep var Ne Cera) vhs easy at eee Stbatetg Pe ars Suk ytce etn rho? aha Ee 44h Pat Pir inta 
Decguniase . ’ sae aie : ry etl Oe a AS Rd gyre see § Ne aN et atane® Ae 
te > ase 8 , aa U r “y yg oe ‘ ae ms gt, 4 Wel A Pir Lang} ere ea rf 7 
J te et “9s yy bee Be ’ eet ts end ! atwejeta’ cd Naa bye y" ¢s ota Cae Se Gee De et ke 
ste aa e seg ats aal ae | aegiagessd at Sy ey a ¢ or Ne Ol aise oc. een POUT Aste 
ue tMee tt ea? ' ' cates bt DK peat ye ag bs, rz Pie he ie] a Meter er eres ty 
%, 8 oe m8 8 seks ; Trt re rig th yt ates nd Seo Tt aesrge ery gle 
J “2 , + et > aa et Abas he i ae sD, 4 i se t's ae ant way Lone 
=s #  # suet i ye < eie . EBhO* wpa ren ee pure ptge ® ae 
r, tee an An “W ara Heaters 1°95 Nerve phere 
+ 2 ‘poh ae an ae i t Aas 10 trg ts 4 oe) Arag 08 aes v4 nee 
seg (stay ee repens iste: 828 420M Fe 
\ > ree he * » 6 ae A Oreos ta 2197 eg tal gt 
ose Cay y au aes , Af lak dha ‘ F4rt et of passes 
, ary Oras ea") ‘ BREEN 9 Tyee ebay tite ee par ee ove 170% oat are 
‘ ComeP reo ‘ GU Eee ee Age eh ees 8 ye Fy trys dade soy nee Hake: a7 9 yeu urte 
Lg ey he hl py fio am ree, tg 4 hd 









Si fitatatel 
tate ate gretyte: mug 

tater cpaesytatetet ane got 

Pe eo ee be ne 

sect geet peak ee eeetats trey tet 

og Chk det dee de 

a aes cf 










Ce) ee Re 2 *yentn pa tapas BEd "Why ea ae t ma ete 
9 * - o 
o's oe koe ay bia.) iad i eae ue 
474 OP i Pear | a bahgte be pe 
























1 e'y4eu ?’ ave 
eg teen on § ter fa bay ata Sete & peu Gee) 
neh MAO, Ga ee a rad ieee 

ta gre TOL OTg 8 aks 


aoa € ite. 4" gignt ats 
* 























: ap 
qe Ween eats 
ar Mut? 
aaa eye? ais 7G oP Lr ie 
‘ rss 






Veh Ovige Sala 
gb es! ty 


oer get, 





fae Dan 





































































































































































































































































pis ¢ 744 és ” . ’ Te 8 ho pathher dat a 2 
. on soy re 73 ‘ ant Ge oe ote ate Saute, mu 2 
oe 4 Le rca fama be rr a ore Ly meee) 4 a ptie. < 
é . “#4 Tou Uber DIME car tat Be eas ew ' 
PC ee tf b Aaa bee ae Sore aap ete it 343- 
es Let he a aaoeee .. Ss 
' a we ae EE Satay nt sm eae) Ci iiytee: Dee re petit ttt tee Favre geaw at ait 
, . . * ‘ . + aint I eh pee 5 
mye Ps oegr hy a ; s¢ ms oe r oat! they wa Pee se itet ee. Hoey eer ele teenies 
‘ an ‘ 2 a.g.0 . i * 
“ ES er | Se ad a 4 “As eee waee AY aoa 2 ealaighe 
é + tA st i seer} 
' mere aer et ba Cher oo ree Y cme anal tt ge 
148 ' apd Pad ‘ t . 
’ Pa etng a pte yah t eel plac it west e ra tire Satta tgs 
Te 5 ba = ta ty at be fie rar Hass 
' ee ee Rywt ee Ps t eM ats nae te ehat ete of . Padite etary seed SF 
Pet me bf ie tos armas tue 1048 3595 yee Adele 2s fing? ‘hs 4 
a6 re pare fou ar) ores hex sugtt a4 rr | Ce gre ies Wu har ee ey oP a tall 
At 24m Wha toe ack, Poa a ¢td jes Ast ay 3 ot pire ge ; CA ha ae » oP yr ahietys oyagmaks 
. be PON ae tere ‘Ter pte ah oe eo rae : Liha et | eheew ites Ds as esy eS 
Sie Orgs go Ueber et anes as cane refye 8 dee, ip , Bre ri eb Pe Lr Mor rd anrurangraruheet ’ 
ie, Wk Heo oe i sh mega y ett thy WV ry detateela fpyatas ere Fy Lad ty ShyZesery 
Peet Sue. acne nan ‘ raat pe eye i set hon esstate uy 2 pow bed 4: jee ait aha 
ati ee Liens Fak Ratio AD "pinta pthae 6m Hipsieie F ‘ Bienes uy LL ae Anta 
. , eet began gees Saottacs Ny ‘ 
auwda « epee the 4 teas 
pihie eg ae er arms Pret? | 
sea + Hag G be gage ac «owt h barrel 
! 2+ a's ae . ' ’ ar pee. 
= ’ . ant Z 
xi ‘ 2 ” Sais re ee I wreath dah gy tty ae 
tos Fp eres os eee chs ‘ ar a at 
he as . ove rine . 
a oye Lt res en » sa » 
Pry + ats Ls hae * t Le 
leet es yeti, oi ec ar (en | 3 aeaad ga fs 
oy Hn ‘ : P * 5 ite 3 ‘a we 
=) se) weg os ae is ev “te oh 
rs PEE SRE > an a atgts * eh? ? em $e Fghotetes 
oe o% tee % tee sNa ¢ i. é i a e ' 
Se tenence gn Oe * Bi ae Ay s aro aU ene states t 
RES eed Heb i AO a REE 
Lae ot 4 op “ge e a 
Pe es “tar tte th? alana 
ie Aes ay 3! 82. Bhar aayisss “ao ew Fae 
*,* 4, bol + heat ‘Gere uz Sie RRL AR ag tr Lite hy ba be Gea 
a ihe ESay sy tng 2° at mais 
. ‘. als yA é De bl 
ae Ee ys: at ' ane os 2 aa 
a” poset? ey rs é tha th | 
«hve > aye teu atc . an) aa : Pa ot? ety tee 
. = Ap yetiag th euigs ane ener ce ‘ if (Meke | 
“455 *> oa wet baat seb cred mah oH h Pi Sor td at aunt = 
. eore ', ty ee ST ee Bt aA * sf 
oe digtin ay era) aaa aS CO tee iat dot al be - Para eect 
a" p> ote @ eee ty 1 Pees wae . . PR tee oh hUee Ee . 1! + 3 2 
® ’ ta fe , s 4 (eo bho ay Tan var pT he Cay. 
' a it, “evptign mipls@als eniecse “ Chae ae" fea sce Sin careeon el 
re ° i Porras rte eee oe eee oy 42,72" vere preter §. 





DUDLEY KNOX LIBRARY . 
NAVAL POSTGRADUATE SCHOOL = §t 
MONTEREY CA 93943-5101 














classified 
irity Classification of this page 


REPORT DOCUMENTATION PAGE 


Report Security Classification: Unclassified 1b Restrictive Markings 


Security Classification Authority 3 Distribution/Availability of Report 
Declassification/Downgrading Schedule Approved for public release; distribution is unlimited. 


~erforming Organization Report Number(s) 5 Monitoring Organization Report Number(s) 


Name of Performing Organization 6b Office Symbol 7a Name of Monitoring Organization 
ival Postgraduate School (if applicable) AA  |Naval Postgraduate School 


Address (city, state, and ZIP code) 7b Address (city, state, and ZIP code) 
onterey CA 93943-5000 Monterey CA 93943-5000 


Name of Funding/Sponsoring Organization 6b Office Symbol 9 Procurement Instrument Identification Number 
(if applicable) 


= (city, state, and ZIP code) 10 Source of Funding Numbers 
Program Element No Work Unit Accession No 


Title (include security classification) TWO - DIMENSIONAL BOUNDARY SURFACES FOR AX] - SYMMETRIC EXTERNAL 
<XANSONIC FLOWS. 


Personal Author(s) Waleed Al-hashel 


a Type of Report 13b Time Covered 14 Date of Report (year, month, day) 15 Page Count 
aster’s Thesis From To MARCH 1993 70 


Supplementary Notation The views expressed in this thesis are those of the author and do not reflect the official policy or position 
the Department of Defense or the U.S. Government. 














Sor ghee 


Cosati Codes 18 Subject Terms (continue on reverse if necessary and identify by block number) 
7 Group. |Subgroup. | TWO - DIMENSIONAL, AXI - SYMMETRIC, EXTERNAL TRANSONIC 
—————— J FLOWS, SMALL PERTURBATION BOUNDARY SURFACES. 


_ J 
Abstract (continue on reverse if necessary and identify by block number) 
vestigation of two-dimensional transonic flows is extended to axisymmetric problems. This is of considerable practical interest, 
r example, with regard to missiles or aircraft engines which approximate much more closely bodies of revolution than two- 
nensional bodies. The main concern with axi-symmetric flows lies not only with the complexity of the governing nonlinear 
rtial differential equation which is mixed of elliptic-hyperbolic type but also with the lack of a general method for accurately 
lving this type of equation. We solve the nonlinear transonic equation using separation of variables technique, which yields two 
nlinear ordinary differential equations. The x-dependence can be integrated numerically, and the solution for the r-dependence 
n be obtained using the expansion method originated by Van Dyke. This works well with only three terms in the expansion. The 
nic solution of these equations is obtained analytically since both equations are simple enough to be integrated for this case (M, 
1.0). The small parameter (1-M,..” ) plays an important role in specifying the shape of the boundary surfaces for external axi- 
mmetric steady flow of interest for design. A Navier-Stokes solver was used to compute the inviscid flow to confirm our results 
the region over the surface where the small perturbations apply. 









Distribution/Availability of Abstract 21 Abstract Security Classification 

€_ unclassified/unlimited __ same as report __ DTIC users | Unclassified 

, Name of Responsible Individual 22b Telephone (include Area Code) 22c Office Symbol 
- Biblarz (408) 656 - 3096 AA/Bi 








FORM 1473,84 MAR 83 APR edition may be used until exhaustedy d security classification of this page 
All other editions are obsolete Unclassified 
ale 


72597930 


Approved for public release; distribution is unlimited. 


TWO - DIMENSIONAL BOUNDARY SURFACES 
FOR 
AXI - SYMMETRIC EXTERNAL TRANSONIC FLOWS 


by 
Waleed Isa Al-hashel 
First Lieutenant, Bahrain Amiri Air Force 


B.S., Northrop University, 1987 


Submitted in partial fulfillment 


of the requirements for the degree of 
MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 


March 1993 oan = 
— 7 


ABSTRACT 


Investigation of two-dimensional transonic flows is 
extended to axi-symmetric problems. This is of considerable 
practical interest, for example, with regard to missiles or 
aircraft engines which approximate much more closely bodies of 
revolution than two-dimensional bodies. The main concern with 
axi-symmetric flows lies not only with the complexity of the 
governing nonlinear partial differential equation which is 
mixed of elliptic-hyperbolic type but also with the lack of a 
general method for accurately solving this type of equation. 
We solve the nonlinear transonic equation using separation 
of variables technique, which yields two nonlinear ordinary 
differential equations. The x-dependence can be integrated 
numerically, and the solution for the r-dependence can be 
obtained using the expansion method originated by Van Dyke. 
This works well with only three terms in the expansion. The 
sonic solution of these equations is obtained analytically 
Since both equations are simple enough to be integrated for 
this case (M,=1.0). The small parameter (1-M,7) plays an 
important role in specifying the shape of the boundary 
Surfaces for external axi-symmetric steady flow of interest 
for design. A Navier-Stokes solver was used to compute the 
inviscid flow to confirm our results in the region over the 


surface where the small perturbations apply. 
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16% INTRODUCTION 


Transonic flow is a classical problem in gas dynamics and 
yet not many exact solutions are available. The main 
Seaecaculty lies, of course, in the nonlinearity of the 
equations and in the boundary conditions. [Ref. 1] 

It is of interest to extend two-dimensional flow solutions 
to axi-symmetric problems, which are more practical for 
designing bodies of revolution such as missiles or aircraft 
engines. The goal is to design a shock-free body or surface so 
we can minimize the penalty of the wave drag that is caused by 
shock waves. The small perturbation technique may be 
implemented for simplifying the governing equations and the 
mathematical procedures for solving the resulting equations in 
nonlinear flows. This research shows a solution to nonlinear 
two-dimensional, external flow, over an afterbody axi- 
symmetric surface. Beginning with the transonic equation for 
an axi-symmetric body in Chapter II and applying separation of 
variables gives two nonlinear ordinary differential equations 
which lead into two exact solutions. One of the nonlinear 
ordinary differential equation can be integrated numerically, 
and other one can be solved using the outer expansion method 
by Van Dyke [Ref. 2]. Chapter III shows the derivation of 


pressure coefficient and local Mach number for different 


surfaces. Supersonic boundary surfaces are obtained 


implicitly in Chapter IV for M=1g905,1.1, and 1.2 .andeene see 


examined by a 3D Navier-Stokes solver, run as an Euler solver 
for this axi-symmetric problem in Computational Fluid Dynamics 
(CFD). In addition, the sonic surface is derived analytically 
in Chapter V. This sonic boundary surface is also examined 


using 3D Navier-Stokes solver in CFD to verify the results. 


dB TRANSONIC EQUATION FOR AXI-SYMMETRIC BODY 
The small perturbation, non-linear, axi-symmetric, 
transonic equation is written as [Ref. 3] 


(1-MZ) 9, + 2 (r9,) = 9,9, (1) 


Nile 


Where the non-dimensional axial, x, and radial, r, velocity 


potentials have been defined as [Ref. 1] 


Z x = 2 a 
= . =e 2 
* , =a (y+1) —, Pp - = 105 (¥ 1) ( ) 


oo 


where M, is free stream Mach number and y 1s ratio of heat 


capacities. 


An exact solution to Equation (1), the transonic equation, 


can be found from 


p (x,r) =@s(xX,7r) + (1-MZ)x (3) 


Which patterned after the two dimensional case given by 


Biblarz [Ref. 1]. Here @* (x,r) is a separable solution for 


the velocity potential. 


The transonic equation in axi-symmetric form, Equation (1), 
may be separated by letting the potential function @ (x,r) 


equal to 


m (X,r) =& (x) 0 (r) + (1-MZ) x (4) 


Substituting the above function into the transonic equation 


results in two ordinary, second order, non-linear differential 





equations 
ae gE ye. d 
ae ERE A— =0 (5) 
q*0 1 a _4 e226 (6) 





where A is the separation constant. 
The solution to the first O.D.E.(5) is obtained by multiplying 


both sides by dé/dx, 
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Thus 


a -( 34¢2 3 (9) 
dx [ 328 +a | 
Rearranging 
dx = ais 
i 10 
(21g +a)? (10) 
2 
and integrating 
d 
AT Ky = : 5 (11) 
(5 46? +a)? 
Z 


where @ and x, are integration constants. 


The solution to the second non-linear O.D.E.(6), 1s obtained 


by an outer expansion method. [Ref. 2] 


C(r) = — + (1-M2) f, (r) + (1-M2)?f,(r) +... (12) 





Where (1-M.) represents a small parameter and the first term 


is the purely sonic solution. 


By taking first and second derivatives and substituting 
them into the original 2nd order, non-linear, O.D.E.(6), the 


solution becomes 


(1-M2) 2 (13) 


(j= 4+|1-Mé|adr ‘v8*2) + atyAzy 2B) 
Ar? 28+8/8 


where ais constant. 


Equations 11 and 13 are solutions to the differential 
equations 5 and 6. Boundary conditions are needed to 


determine the constants a, C,, a, and A. 
The constants C,, a, and A can be shown to be related by the 


expression [Ref.1] 


C 1.2426 : 
ee = 1.08x107 (14) 
A} 40-7574 
and 
Cae ee ood by ed re (15) 


The positive sign above is for M, = 1.0 and the negative sign 


is. tor Mor aoe 


Weonmiunsertingebqn.(15) into Eqn.(11), we have 


KaX,= de 5 
3 2 2 5 
2 AE? 4+ C, |1— Me |2-7579)3 
Factoring out C,, we obtain 
es A 
a= Cine § 5 
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3 Baa m2 [2-757 
eee 


Introducing new variables 





F (V8 +2) =r VAS) Xr a 


Equations (17) and (13) become 


= wat 
x- x5" C,°(38) ae 


(E24 [a m2 |2-7574]3 


(16) 


(17) 


(18) 


(19) 


(20) 


2,2 
¢= Lak) OE a +|a-me |e v2) pve (21) 
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Finally, defining two new variables 


i il 
POG © se (32)? (22) 
2 
ae ae aan 
(a jh) Dem 
Equations 20 and 21 become 
x= ee ie 
(aes Le Ee 2 
where Xp=0. 
(1-M2)* 
‘Co A+ |WeMe | 7 08° 2) 5 eee (25) 
me 50.635 


Equation (24) will be numerically integrated and plotted in 


Fig. (1) as € (%) vs. & , and Equation (25) will be evaluated 


and plotted in Fig.(2) as C(f)ivs. £f for  Mj= 1905 eee 


jee These results will be used later to obtain transonic 


surfaces with their corresponding C, and local Mach number 


peerates. Im Big. (2), Equation (25) is plotted up to its 
minimum value with increasing — Thereafter, a constant 


value is patched in to the right of the + mark. This is done 
because of the need to keep € constant reflecting the 


necessary behavior of the function for range fF [Ref.1]. 


JET GIES PRESSURE COEFFICIENT AND LOCAL MACH NUMBER 


The pressure coefficient for a flow with small 


perturbations is given by [Ref. 4] 


2 2 
2Uy . (y_yzy (Bx) . (Bz 
U oe Op. 


oO 


(26) 











The linearized pressure approximation for axi-symmetric flow 


is 


Cc. = = (27) 





Recall the axial velocity potential, Equation (2), Chapter II, 


yp, = Me (y+1) = (28) 


oO 


thus 


a (29) 
U, M2 (y+1) 
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substitute Equation (29) into Equation (27), yields 


-? ~, 
o = — (30) 
ey te) 


We need to take the derivative of the potential function with 


respect to x in equation (4). 


n= @ SOS a oe 


cee 
aE (31) 


Rewriting equation (9) with the constant a inserted becomes 





dé _[ 32 pg _2) |2-7574]3 (32) 
ie ee +C, | (1-Mé) | 
Recall Equation (18) 
1 i 
= i ( 3 i (33) 
BO. 


Mmoieerting Eqn.(33) into Eqn. (32) and factoring C, out, yields 
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Recall Equation (23) and rearrange 


5 (a Pyc-stec 
ca 


Rewriting Equation (31) the potential function as 


(aye - 


3 i 
Q, = a a 6 ee | 7+] (1-m2) fe ave) 


Thus 
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i 
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Recall Equation (14) Chapter II 


C 1.2426 
a = 1.08x10~ 


then finally Equation (37) becomes 


P= 022208 Cire | (las) Tae 
Rewriting the coefficient of pressure 
Pp 


Mz (y +1) 


Vz 


& -_=2__|.22082(82+ | (1-2) feats hue) 


(35) 


(36) 


(37) 


(38) 


(335) 


(40) 


iMempGevylous expression) for) Cy will be used further mel 


determining the local Mach number. 
The perturbation velocities are assumed to be vary small 


compared to the undisturbed uniform velocity JU,. [Ref. 5] 


Hence at a point near the body, the velocity vector Wh eens 


2-D flow is given by 


aan (41) 


We can obtain the local Mach number M in terms of the 


perturbation velocities and the local speed of sound c¢c as 


Jake ee (42) 
Ga C2 


M2 = 


Neglecting the higher order ratios of the free stream to 
perturbation velocities, since they are considerably smaller 


than unity, we have 





2 2u, 
7 s + U. | (43) 


i 


The ratio of the local speed of sound to the undisturbed 


uniform speed of sound becomes 





ee ie 
or 
c2 _ T,/T. 
Pe T,/T (45) 


Where T, is the stagnation temperature. 


We can write temperature ratios as functions of Mach numbers 


as 





nN 
NO 


(46) 


Aife 


eee (Y-1) ye 
Z 


Substituting Equation (45) into Equation (46) and rearrange, 


we have 


an Z 


Viz ieee 
[a+ tae) 


1+ fa « (Y=2) oe 
M2 U. (47) 
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from Equation (27) 





48 
u (48) 


Rearranging and solving for local Mach number 


2 
y2 = Ma (1-Cp) 


1+ 71 
2 





: (49) 
Me. Cy, 


This equation will be used to show the local Mach number for 


Seren SUuprace at M. = 1.05, 1.1 and 1.2 [Ref. 6]. 


aS 


IV. SUPERSONIC BOUNDARY SURFACES 


The boundary conditions which need to be applied require 
that the gradient of @ vanish far ahead of the body and that 
the flow be tangential to the surface [ Ref. 7]. In terms of 


perturbation velocities the boundary conditions become 


u 
=e (50) 


we 


dr 
ce 


| surface 


Recall the modified velocity perturbation potential Eqn. (2) 


u 
p, = Me (y+1)— (51) 
Us 
then, we have 
u, Y-, 
aL <2 ae (52) 
U, MS (y+1) 


by substituting Equation (52) into Equation (50), we obtain 
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By differentiating Equation (4) with respect to r gives 


ac 
= = 54 
ee dr (54) 


Recalling Eqn.(18) and substituting Eqn. (21) into Eqn. (23), 


and taking the derivative with respect to fF, the above 


equation becomes 


dv 
= 0.0848 —— (55) 
Pr g ain 


Substituting Equation (55) and (13) into equation (53), we 
have 
eS pS =o 2 go84 1M | £2 -9284 


ax Mz (y+1) de (56) 


mde Sia (ll Ms): o>! 


It can be noticed here that € and € are given implicitly as 


function of X and fF respectively. 


From Eqns.(19) and (22), we can rearrange Equation (56) as 


-2 Rs, 
Bee eo end z = + 2.8284 Ja-m2 | #1.8284 


Mz (y +1) a (57) 


8.18 


Or woe Ma)? 7 oS! 
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Further arrangment of the previous equation yields 


df _ 
20 8264 es pi-0284 _ 9 91512 (1 =i) ee 
Ye (58) 
= =2 
—— E dx 
Mz (y¥+1) 


Now, we can compute the exact boundary surfaces out of the 


equation above. Starting with the left hand side of the 
equation 
ae ie 
L(f) -—_____- (59) 


[8-2 8284 |1-M2|F4-8284-9 .1512 (1-M2) aad 


By defining equation (59) we can plot Z (r) vs. ¥F in) Fage 


(2a sE OG ia) 4a10 5 oe te cool ola. 
The minimum values of fF, for different Mach numbers can be 
expressed as fF, = F,(M,) which has been found to be [Ref. 1] 


oO 


2 ere 


O Ja —m2|°-2072 (60) 


als: 


Fig. (4) represents fF, vs. M, and shows a symmetry close to 
fe-1.0 where fi=o. So for our purpose to represent the 
condition sufficiently close to M, = 1.0, we will choose a 
finite value of fF, (perhaps as 4.0). 


Defining the right hand side of Equation (58) as 


= =2Z ‘ 
K() oad F ax (61) 
Ms (y¥+1) 


taking the derivative of Equation (24) and substitute it into 


Equation (61), yields 


g 

_ -2 

rire se oe | Edt - (62) 
Mes victels) [E24 |1-m2|2-7574] 3 


integrating by parts, we obtain 


Z 2 
K(2) = =22 REY | @2— [ane p-7576)3 [a -z pear (63) 
Carers 


on? 


the equation above will allow us to plot K(X)vs.X in 


Bree. 5) for M, = 1.05, 1.1, and 1.2. Where K(X) < 0 


femeiMw= 1.0 and K(X)> 0 for M, <1.0. 


abe 


Therefore, integrating both sides of equation (58) becomes 


AY 


[= -2 .8284 |a -m2 | #2-8284 _ 1512 (1M 
sl x (64) 
x 


MZ ( yas 
a Yom 


Equation(64) is numerically integrated and which then allows 
us to determine the boundary surfaces in dimensionalize and 
non-dimensionalize (normalize) form as shown in Figs. 6 and 7 


for M, = 1.05, 1.1, and 1.2. Now, we can determine Co from 


Eqn. (40) and local Mach number from Eqn. (49) for the three 
surfaces obtained at different Mach numbers. 

In Figs. 8-13, the supersonic boundary surfaces for 
M,=1.05,1.10, and 1.20 are plotted with their pressure 
coefficients and local Mach profiles in dimensionalized and 
non-dimensionalized (normalized) form. These results will be 
verified later using a 3-D Navier-Stokes solver in 


Computational Fluid Dynamics (CFD). 
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SONIC BOUNDARY SURFACE 


The resulting equations for 


can be derived in explicit form 
(65) 


The sonic flow (M,=1.0) 
(24) and (25). 

z3 
(66) 


from Equations 
sonic flow are 


and 

4 

Ee 
(67) 


Rearranging Equation (65) 
is integrable, we obtain 


Since Equation (67) 


Or 
21 


(68) 


(69) 


We can write Equations (11) and (13) in sonic form as 





d 
x = — = (70) 
3 2\3 
(7 2) 
and 
(<— (71) 


integrating Equation (70), we have 


=. eal 3 
= = 72 
leg, Ts A xX (72) 
From Equation (4), we can write the sonic perturbation 
potential as 

p(x, ©) = — Coe) (73) 


Substituting Equations (71) and (72) into Equation (73) yieree 


IE Ge 
,r) === 74 
p(x, TF) er (74) 
Equation (57) can be written for M,=1.0 as 
ar. = 0.2608 & (75) 
Ox (y +i) #3 


Ze 


Rearranging Equation (75), we have 


£3 dF =-0.1087 E dz& (76) 


Therefore, integrating both sides of Equation (76), we obtain 
f 

[sar =-o.1007 [ E dx (77) 
fe 


Inserting Equation (69) into Equation (77) yields 


74 -F2=-0.0040 X4 (78) 


where f,=4.0 has been chosen as an asymptotic value to 
represent the condition sufficiently well at M,=1.0. Equation 


(78) will allow us to determine sonic boundary surface. 


The pressure coefficient can be written from Eqn. (40) as 
_ -0.4416 975 
C, === ( £? (79) 


in case of M,=1.0. 


Ze 


The local Mach number can be written from Equation (49) for 


Ehe sonic sSsurEace as 


M2 = (1 - C,) (80) 
b+ >| 
Ss Pp 
The sonic surface is represented dimensionally in Fig. (14) 


with the pressure coefficient and local Mach number profiles, 


and normalized in Fig. (15) with 7,=4.0 for (M,=1.0). The 


sonic boundary surfaces, in both forms, are plotted with the 


Supersonic boundary surfaces in Figs. 16 and 17. 
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VI. COMPUTATIONAL FLUID DYNAMICS (CFD) 


The history of Computational Fluid Dynamics (CFD) is 
closely tied to the rapid advances in digital computer which 
has a great impact on problems of design in modern engineering 
practice. These problems can now be solved at very little cost 
in a few seconds of computer time which would have taken years 
to work out with the computational methods and computers 
available twenty years ago. [Ref. 8] 

The CFD will be used to compute the axi-symmetric flow 
over the boundary surfaces obtained by the small perturbation 
method. 

The computer program GRAPE [Ref.9], an acronym derived 
from Grids about Airfoils using Poisson’s Equation, has been 
written by R. Sorenson at Ames Research Center to generate 
two-dimensional finite difference grids about airfoils and 
other shapes by the use of the Poisson differential equation. 
Outer and inner boundaries are specified for the C-type grid, 
where the surface of the body is treated as the inner 
boundary. Important characteristics in grid generation are 
the ability to specify the spacing between mesh points at the 
boundary, in the direction normal to the boundary, and the 
control of the angles with which mesh lines intersect the 


boundaries which is known as orthogonality. 


Z5 


The grid size for the sonic surface (M,=1.0) is 107 x 60 and 
for the supersonic surfaces (M.=1.1, 1.2) is 115 x 60a 
shown in Fig.(18). Since the surfaces are symmetrical, half of 
the grid or eventually the lower surface was rotated for 11 
planes plus two more for symmetry to generate an axi- 
symmetric after body surface. 

The OVERFLOW program [Ref.10] developed by Ames Research 
Center which uses 3-D Navier-Stokes and Euler solver for 
viscous/ inviscid flow was applied. The results are shown in 
Figs. 19-21 for three different boundary surfaces at M, = 1.20; 
1.1, and 1.2 with 1500 iterations. It was found that the 
density residuals decreased by more than one order of 
magnitude over 1500 iterations. After 3000 iterations, the 
solution converged by two orders of magnitude. 

In case of the sonic surface Fig.(19) with M,=1.0, the 
Mach lines contour shows the shock is forming downstream at a 
local Mach M=1.65 at which the flow becomes’ subsonic 
downstream of the shock. This result complies with the small 
perturbation transonic solution obtained earlier in Fig.(15). 
In other words, the flow is shock free over the sonic surface 
ieee Med Ol 

For the supersonic boundary surfaces M.=1.10, and 1.20, we 
can see that the flow starts at M=1.10 and M=1.20 respectively 
and follows the afterbody surface until it shocks. The shock 
is formed at local Mach M=1.75 and 1.85 respectively as shown 


in Figs. 20e@and/21- 
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On the other hand, these results agree with the small 
perturbation solution over the transonic range as plotted in 
Figs. 11 and 13 where the local Mach is represented by the 
dash line and it is increasing towards the afterbody surface 
until it shocks. This indicates that the flow is shockless 
over these supersonic boundary surfaces in the transonic 
regime. 

The velocity vectors are plotted in Fig.(22) for M,=1.0 
which shows that the flow is following the afterbody surface 
until the shock, and then flow separation takes place due to 
the steep pressure rise and the introduction of numerical 
viscosity into the flowfield at this location. 

These boundary surfaces obtained are of interest for 
design of body of revolution such as missiles afterbody or 
aircraft engines, in which a patching technique may be applied 
to these surfaces to get the best shockless surface in the 


transonic range. 


Zh 


Vit. CONCLUSION AND RECOMMENDATION 


We have shown in this research a solution to the non- 
linear transonic small perturbation equation by implementing 
the separation of variables technique. The x-dependence was 
integrated numerically and the r-dependence was solved by an 
outer expansion method. It was found that the parameter (1- 
M7) has strong effect in specifying the shape of the boundary 
surfaces of interest for design. A 3-D viscous/ inviscid 
Navier-Stokes/ Euler solver in Computational Fluid Dynamics 
(CFD) confirmed our results obtained from the small 
perturbation technique for the boundary surfaces at M,=1.0, 
1.10, and 1.20 . In other words, we achieved our goal of 
designing shockless surface in the transonic range. 

This research may be extended using a patching technique 
to obtain the best shockless surface in the transonic flows. 
Patching criteria will have to be developed based on the 
mathematics and the flow constraints. Also, connection between 
small perturbation solution and CFD might eventually be done 
internally in the computer. The small perturbation is the most 
efficient method to define or to design axi-symmetric 
afterbody transonic surfaces, however CFD can be used to 


predict the performance of these aerodynamic surfaces. 
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APPENDIX B 
COMPUTER PROGRAMS 


KEKE KKKKKEKKKEK KEK KEK KEKE KKK KKKKKKEKEKKKEKKEKKEKKKEKKKKKKKKKKKEKKKKKK 


* 
* 
* 
* 
* 
* 
* 


100 


This program is designed to calculate KSI(X), 
K(X), and X using numerical integration based 


on trapezoidal rule to solve Eqn. (24) and 
plotting the output as shown in Figs.1 & 3 


KKEKKKKKKKKKEKEKEKKEKKEKKKKKEKKKKKKEKKKKKKKKKKKKKKKKKKKKKKEKKRKKKRKKEK 


PROGRAM KSI 

REAL M(3),%(0:401,3),A,A1,B(401, 3) 

+ »KSI,P, FUNC, H, DD, O,L,K(401, 3) 
INTEGER YY 

OPEN (UNIT=9, FILE=’ KSI’ , STATUS=’ UNKNOWN’ ) 


PRINT *,’ENTER LOWER & UPPER BOUND, # OF INTERVALS, 


+#OF DATA SET’ 
READ *, A, Al, N, Ni 


DO: TO =i 33 

PRINT *,’ENTER MACH NO. ' 
READ *, M(I) 

M(I) = 0.932705 


IF (M(1) . Lelio ann 
Jeo ike 

ELSE 

P=1. 

ENDIF 


DO 20 J=1,N1 

B(d, 1) =o0*f09 il 

H = (B(J,1)-A)/N 

AREA = 0. 
K(J,1)=(-0.0102/7 (M(1) **2)) * ( (ABS ( (Big ieee 
+( (ABS (1-M(1)*220pee nd. 7574) 26 8 ee 

+ (ABS (1-M(I) **2) )**1.1716) 


DO 30 L=1 , N 
KSI = A + (L-0.5)*H 
DD = (KSI**2+P* (ABS (1-M(I) **2) ) **1.7574) 
IF (DD.GT.0.) THEN 
FUNG = 12/7 (Bal 3) 
0) Seale 
ELSE 
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* 
* 
* 
* 
* 
* 
* 


FUNC = 1./((-1*DD) **(1./3.)) 
Q= -l. 
ENDIF 

AREA = AREA + H*FUNC*Q 

X(J,I) = AREA 


30 CONTINUE 
PAY CONTINUE 
186, CONTINUE 


DO 40 J=1,N1 

Peeniese, (B(I, Ly xi ,1),K Gel), l=1,3) 
50 HORM ix, 3(LF6.S 7250.5, 1xX)) 
40 CONTINUE 


200 PRINT*,’TRY AGAIN ? ENTER 1 FOR YES,2 FOR DATA 
+FILE,OTHERS FOR NO’ 


READ *, YY 

fee (YY.BO.1) THEN 
GO TO 100 

ELSE 


IF (YY.EQ.2) THEN 
DO 0 «d= 11 
WRITE (9,50) (B(J,1I),X(J,1I),K(J,I),I=1,3) 
AL ALG, CONTINUE 
ES) TS). Bele 
ENDIF 
ENDIF 
END 
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KKH KK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKKEKKKKEKKKKEKEK 


* * 
* This program is written to calculate ZETA(r) and * 
* rousing Egqn.(25) and plotting the output oto * 
* * 


KKEKKKKEKKEKKEKEKKKEKKEKEEKEEKEKEKKKKEKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 


program zeta 
real M(3),zeta(0:14,3),r(14) 
open (unit=9,file='’zz’',status=’ unknown’ ) 


Gem i.0) L=4e43 

c M(I)=0.7+1I*0.1 
Print *," enter macheno 
read *, M(T) 


ado 20) JV=i,14 
zeta(0,1I)=1000 
print *,’enter r values’ 
read *, r(J) 

C me) = Oe sa) 
zeta(J,1I)=(4/ (r(J) **2))+ ( 
+(r(J)**2.8284) )+((1- (M(I)* 
if (zeta(J,I) .GE.zeta(J-1, 

20 continue 

ale, continue 


abs (1- (M(I) **2) )* 
*2))**2)*% (r (J) **7.657) /50.63 
I)) zeta(J,1i)=zeta(J-1,1) 


pDrlat 30). eM 
30 format (11ix,SF10.4) 


ado 40 Jai 14 


c Print. 50, -2(o)-(2erat hae 
write (9,50) (eid) »Zetatu )). i— ies) 
50 format. (4% .3:(2FP10. 4.29) 
40 continue 
end 
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* * 
* This program 1S written to calculate 2Z(r) and r * 
* using Eqn. (59) and showing the output in Fig.4 . * 


* * 


KEKE KKEKKEKEKEKEKEKEKKEKEKKKEKEKKEKKKEKKEKEKEKEKKKEKEKKEKEKKEKKEKEKKKKKKKKKKKKKKK 


PROGRAM ZR 
REAL M(4),R(1000) ,Z(0:1000, 4) 
OPEN (UNIT=99, FILE=’ ZR1’ , STATUS=’ UNKNOWN’ ) 


DO 10 I=1,4 

€ Mier) = 0.7 -+ 1*0.1 
PRINT *, ‘ENTER MACH NO. ’ 
READ *, M(T) 


Ber20 J=1,1000 
miO.1) = 0.0 
R(J) = 0.001*J 
Z, ae Sle) = (Ree tees X30 see ( 8 
(2.8284* (R(J) **4.8284) *ABS (1- (M(I) **2)))- 
+(0.1512* (R(J) **9.657) * ((1-M(I) **2) **2))) 
mamaigd,l).LE.0.0) Z(J,1I)=Z(J-1,1) 
20 CONTINUE 
10 CONTINUE 


DO 40 J=1,1000 


Cc pean 50, RJ) , (Z2(J,1),1=1,4) 
Pane (99750) R(J), (Z2Z(0,1), I=1, 4) 
50 FORMAT (4X,5F10.6) 
40 CONTINUE 
END 


oo 


KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KK KEKEKEKEKEKEKKKKEKEKKEKEKKKEKEKKKKEKKKKKKKS 


* * 
* This program iS written to determine Rmin for * 
* transonic Mach numbers using Eqn. (60) and the * 
* output is shown in Fig.5 * 
* * 
* * 


KK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KK KKK KEK KEKE KKKKKKKEKKKKKKKKKKKSE 


PROGRAM RMIN 
REAL M,R 
OPEN (UNIT=88, FILE=’ RMIN’ ,STATUS=’ unknown’ ) 


DO! 200 M = 0-6 12240- 0000 
ITF(M.LT.1.0) THEN 
R=1.2074/ ( (1-M**2) **0.2071) 
ELSE 
TE (MeGTr. 1.00) Dean 
R=1.2074/ ( (M**2-1) **0.2071) 


ELSE 
R = 7.0 
ENDIF 
ENDIF 
c PRINT 80, M,R 
WREERE (887,80) M,_R 
80 BORMAT (ih. 5, 7h lr oeG) 
100 CONTINUE 
END 
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KKEKKKKKEKEKKEKEKKEKKEKKKEKKEKEKKEKEKKEKKEKEKKKEKKEKKEEKEKEKKKEKEKKEKEKKEKKEKEKKEKKEKKEK 


* 
* 
* 
* 


This program is designed to numerically integrate 
Z(r) ,Eqn.(64), based on trapezoidal rule 


* 
* 
* 
* 


KKEKKEKKEKEKKKEKEKEKKKEKEKKKEKEKKKKEKEKKKEKEKKKEKKEKKKEKEKEKKKEKEKEKKKKEKKKKKKKKKKKEKEKK 


100 


PROGRAM IZR 


REAL M(3),R(401,3),A,A1,1Z2(0:401,3) ,H, FUNC,N,N1,RO 


MNTEGER YY 


OPEN (UNIT=55, FILE='IZR’ ,STATUS=’ UNKNOWN’ ) 


DO 10 T=1,3 


M(I) = 0.7 + I*Q.1 
PRINT *, “ENTER MACH NO. 


READ *, M(I) 


PRINT *, ‘DETERMINE THE VALUES OF MACH #’ ,M(T) 
PRINT *, ’ENTER Rmin 


READ *, A 


Pret *, ‘ENTER THE UPPER LIMIT OF R ,Al’ 


READ *, Al 


PRINT *, ‘’ENTER # OF INTERVALS 


READ *,N 


PRINT *, ‘HOW MANY DATA SET DERIVED ? N11’ 


READ *, Nil 


PRINT *, ’INPUT Ro FOR THIS MACH # 


READ *, RO 


Oreo J=1,N1 
Rigen) = J*A1/N1 
H = (R(J,1I)-A)/N 
AREA = 0. 


DO 30 L=1,N 
RR = A+(L-0.5) *H 


FUNC= (RR**3) /(8- (2.8284* (RR**4. 8284) *ABS (1- (M(I) **2) )) 
+ -(0.1512* (RR**9.657) *((1-M(I) **2.)**2.))) 

AREA = AREA + H*FUNC 

IZ(J,I) = Ro - AREA 


CONTINUE 
CONTINUE 
CONTINUE 


DO 40 J=1,N1 


eee (kh io,0),L2(0,1),il=1,3) 


a 


S/ 


a 


Re. 


3,18, 
40 


WRITE(55,50) (Rid ,2),12(9 723 ae 
PORMAT (135) 3 P6 eon Bo cos 
CONTINUE 


PRINT *,’'TRY AGAIN ? ENTER 1 FOR #YES#, OTHERS FOR #NO#’ 
READ *7 12 


fe (YY.EO. 1) Goro sae 
END 
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KEKKKKEKEKEKKEKEKEKKKKEKEKKKEKKEKEEKEK KEK KKK EKKKEKEKEKKEKEKEKEKKEKKEKKKKEKKEKKKKEKKEKKKK 


* * 
* This program is using numerical integration to * 
* determine the supersonic boundary surfaces , a 
* Eqn. (64), and calculating pressure coefficient * 
* and local Mach number profiles, using Eqns. (40 * 
* & 49). The output is shown in Figs. (6-13) * 
* * 
* * 


KEKKKKEKEKEKKKEKEKKKKEKEKEKKEKEKEKEKKEEKEKEKEKEKEKKEKEKEKKKEEKKKEKEEKEKKEEEEKEKEKKEKEKEEKEKE 


PROGRAM CPP 
REAL X(401,9),Y(401,6), XR(401) ,YR(401) ,M,MX, XYMAT 
OVCEGER £,0,K, lL, IWRITE 


OPEN (UNIT=52, FILE=’cp15.dat’ , STATUS=’ UNKNOWN’ ) 
OPEN (UNIT=53, FILE=’cp11.dat’ , STATUS=' UNKNOWN’ ) 
OPEN (UNIT=54, FILE=’cp12.dat’ , STATUS='’ UNKNOWN’ ) 
OPEN (UNIT=77, FILE=’cpp.dat’ ,STATUS=’ UNKNOWN’ ) 
OPEN (UNIT=78, FILE='m105.dat’ , STATUS=’ UNKNOWN’ ) 
OPEN (UNIT=79, FILE=’tirel.dat’ , STATUS='’ UNKNOWN’ ) 
OPEN (UNIT=80, FILE=’tire2.dat’ , STATUS=’ UNKNOWN’ ) 


100 OPEN (UNIT=1, FILE=’ KSI.DAT’ , STATUS=’ OLD’ ) 
OPEN (UNIT=2, FILE=’ 1ZR.DAT’ , STATUS=‘ OLD’ ) 
REWIND (UNIT=1) 

REWIND (UNIT=2 ) 


IWRITE =78 
MACP=52 
ReaD eI *, END=50) ((X(L,J) ,J=1,9) ,I=1,401) 
58 READ (2,*,END=60) ((Y(K,L) ,L=1,6) ,K=1,401) 
60 Wee Ji—2,6,2 


PRINT *,’ Ro is rmin ; EPSILON is the accuracy ' 
PRINT *,’ENTER MACH NO., Ro , EPSILON’ 

READ *, M,RO,EPS 

WRITE (77,*) M,RO,EPS 

AM = 1-M**2 

L = J/2*3 


me 30 I a4 OF: 
DOS ZO0 K 1,401 
Reel, ) ol .0.001) THEN 
Conroe 10 
ELSE 
XYMAT = ABS(Y(1I,JdJ) -ABS (X(K,L))) 
IF (XYMAT.GT.EPS) THEN 
GO TO 20 
ELSE 


Hil 
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XR K) = Xe Re 
YR (1) ey ( fwd) Re 
ZETA = (4./(Y(I,J-1)**2)) + ABS (XM) *Y(1,J-1) **2yGeee 


+ XM**2a0(Y (1, 0-1) #27 65¢ny Somae 

Cp = (-2/(M**2*2.4))*(0.2208*ZETA* (X(K,L-2)**2 + 
(ABS (XM) ) **1.7574) **(1./3.) +XM) 

MX = SORT (ABS ( (M**2* (1-Cp) /(1+0.2*M**2*Cp) ))) 

PRINT 90, X(K,L-1),Y(I,d-1),Cp,MX 

WRITE(77,40) X(K,L-1),¥(1,0- a) Goer 

WRITE (MACP,40) XR(K) ,YR(I) ,Cp ,MX 

WRITE (IWRITE,70) XR(K) , YR(I) 

ele) avo) Se 

ENDIF 
ENDIF 


CONTINUE 
CONTINUE 
IWRITE=IWRITE+1 
MACP =MACP+1 
CONTINUE 


FORMAT (1X,4(F15.7,2X) ) 
PORMAT (1x) 2Pi7 7.9) 
FORMAT (1X,4(F15.7,2X) ) 


PRINT *, TRY AGAIN ? 1 FOR YES!,2 FOR NOW 


READ *, IY 
HEY BO.) (CO vo 100 
END 


60 





KKEKKEKKEK KKK KEKE KKK KEKE KEKE KKK KEKE KEKE KEKE KKK KEKE KK KKK KE KKK KEKE KKKEKEKEKEKEK 


* * 
* This program is written to determine the sonic * 
* boundary surface M = 1.0, ro = 4.0 using Eqn. (78), * 
* and calculating pressure coefficient and local * 
* Mach profile using Eqns.(79 & 80). The output is * 
* shown in Figs. 16 & 17 . Some output files are used * 
* shel (Lye * 
* * 
* * 


KEKKEKKKEKKE KEKE KEKE KKK KEKE KK KEKE KKK KEKE KEKE KEKE KK KEKE KEKE KEKEKEKEKEKKEKEKEKEKEK 


program sonic 

real M(0:160) ,r(0:160) ,x(0:160) , zeta(0:160) ,ksi(0:160), 
+ Muro) a0. iGO), €o(0:160), RR(0:160),YY(0:160), 
+ xx (0:160) ,xxx(0:160) 


open (unit=11, file=’sonici.dat’,status=’unknown’ ) 
open (unit=22, file=’sonic2.dat’,status=’ unknown’ ) 
open (unit=33,file='’sonic3.dat’,status=’unknown’ ) 
open (unit=66,file='’sonic4.dat’,status=’unknown’ ) 


@ ZoMUroumGdat ==> x , ©, Cp , M 
e sonic2.dat => x , r (horintal) 
C sonic3.dat => x , r (vertical) 
Cc POMmied. dat => <x/ro , £/ro , Cp , M 
Goeto £f = 0, 159 
pe(~.) = I*0.1 
ro = 4.0 
r(I) = (abs(ro**4. - 0.004* (x(I)**4.)))**0.25 
RR(I) = r(I)/ro 
YY(I) = - 1.0* RR(T) 
pex(l) = x(1I)/ro 
xxx(I) = 3.975 - xx(TI) 
zeta(I) = 4./(r(1I)**2.) 
eset) = (x(~)**3.)/27.0 
mud) = ri{l)**3. 
k(I) = - 0.00101* (x(I)**4.) 
Co(I) = - 0.1840 * zeta(I) * (ksi (I) **0.6667) 
eee eos (teat) / (1 + 0.2*Cp(T))))**0.50 
S: fines *, xx (I) ,RR(I) ,Cpo(I) ,M(T) 
sae continue 


do 20 I =0, 159 

Wome i), 50) xx(I) , rr({I) , Cp(I) , M(I) 
prette (22,60) xxx(T) , RR(I) , YY(I) 
meite(33,70) xxx(1),YY(I) 


oa 


80 


30 


write(66,50) x( 1) aa eo em 
format (2%,1£8.4)72x%) 384 
format (2x, iui gea2, Zee eee) 
format (10x, 2842. 6 Se hee 
continue 


write(34, 80) (xxx 
write (34,80) (xxx 
write (34,80) (yy 
write (34,80) ( 
EOrMat(s (1x, £2 


(Be 10) aS) ) 
(1) ,1=156,0,-3) 
(1) ,1=0,159,3) 
are = 18s 156,0,-3) 
Ze ‘ae a 


do 30 I =158, 0 ,-1 
write(33,70) xxx(I1),RR(I) 
continue 

end 
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